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ABSTRACT 

We have modelled ultracompact Hll regions (UCHuR) in terms of steady subsonic 
ionized flows in a clumpy medium. Mass loss from neutral clumps allows the regions 
to be long-lived. We examine the form of global flows for different dependences of 
the volume mass injection rate, q, on radius and Mach number, and describe the 
solutions in detail. We find that three observed UCHuR morphologies are reproduced 
with these models. Mach number independent flows that include a radial variation 
can give centre- brightened core-halo morphologies. Mach number dependent flows 
reproduce naturally the uniform UCHuR morphology. In a hybrid model, including 
subsonic and supersonic flows, we allow a supersonic wind to shock in the ionized 
region. The ionized subsonic gas has a high density and so dominates the emission. 
The shell produced has a velocity structure very different from that of fully supersonic 
models. Several morphologies of spherical UCHuR can be understood in terms of these 
various models; however, kinematic data are crucial as a discriminant between them. 
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1 INTRODUCTION 

The importance of massive stars in the Galaxy is hard to 
overstate. Over their relatively brief lifetimes they inject 
large amounts of energy and momentum into the interstellar 
medium through their UV radiation fields and high- velocity 
winds. Eventually they explode as supernovae, adding heavy 
elements, and driving Shockwaves into their local environ- 
ment. To understand the global properties of our Galaxy a 
detailed knowledge of the evolution of massive stars is es- 
sential. Because of their powerful winds and radiation fields, 
OB stars soon disrupt their natal environment, making it 
difficult to observe the conditions in which they formed. 



Howeve r, ultracompact Hll regions (UCHuR) arc promis- 
ing objects of study as they contain relatively young Ub 
stars that have not yet dispersed the original cocoon of gas 
from which they formed. 

UCHuR are embedded deep within molecular clouds 



((n1)L Si 10 7 cm 6 pc). L is the distance along the line of 
sight in parsecs and n c is the number density of electrons. 

Wood & Churchwell fll989| ) proposed that UCHuR mor- 
phologies could be described by a few types: cometary (~ 20 
per cent), core-halo (~ 16 per cent), shell (~ 4 per cent), 
irregular or multiply peaked (~ 17 per cent) and spheri- 
cal/unresolved (~ 43 per cent). UCHuR vary greatly in ap- 
pearance, and the assignment of a particular morphological 
type may, in many cases, be rather subjective. In particular, 
surveys take snapshots of individual regions at only one or 
two different spatial scales, and interferometers such as the 
VLA can miss large-scale structures. Recent long-d uration, 
multiple VLA configuration obs ervations of Sgr B2 ( Glaumc 



et al. 1995 



De Pree et al. 1995) dramatically illustrate the 



and an; obscured by dust so that they can only be ob- 



complicated interactions between stars and their local envi- 
ronment. 

Most attention has been directed towards t he cometary 
UCHuR. The steady-state bow shock model ( Van Burcn 



et al. 199C 



Mac Low et al. 1991 



Van Buren & Mac 



served in the radio and far-infrared (FIR) wavebands. Low 1992| ) has been questioned recently. For example 



At these wavelengths, they are amongst the brightest 
compact sources in the Galaxy. The typical characteris- 
tics of UCHuR, as reviewed by C hurchwell (199C) and 
Kurtz, Churchwell fc Wood (}l994|, hereafter KCW), arc 



two 

well-known cometary UCHuR (including the prototypical 
cometary UCHuR, G34.3 + 0.2C) have been shown to have 
a tri-limbed tail structure with a l arge velocity gradient per- 



that U 2HllR have small diameters 10 17 cm , are dense 
({rie) 1 / 10 5 cm~' :i ) and have high emission measures 



pend i cular to the head- tail axis (Gaume, Fey & Claussen 
1994; Gaume et al. 1995). These observations are incompat- 
ible with the bow shock model. Moreover, the bow shock 
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model cannot reproduce the many arc-like UCHuR seen in 
the new image s. We model such arc-like UCHuR in a forth- 
coming paper ( Williams, Dyson fc Redman 1996 , Paper III) . 

Dyson (L994) suggested that the clumpy nature of 
molecular clouds could account for the relatively long in- 
ferred lifetimes of UCHuR. The clumps act as localized 
sources of mass which is added slowly to the flow by pho- 
toionization and/or hydrodynamic ablation. This continu- 
ous mass injection leads to a recombination front bounded 
UCHuR that does not expand quickly, and so does not lead 
to the lifetime problems encountered if UCHuR are mod- 



elled a s 'classical' Hn regi ons. In a previous paper (Dyson 
Williams & Redman 1995, hereafter Paper I) , we examined 



one of the simple models outlined by Dyson ( 1994 ). We cal- 
culated line profiles and emission measures for a supersonic 
wind driven UCHuR. In this model, the ionized flow remains 
supersonic through to a recombination front. This type of 
solution reproduces the shell morphology of some UCHuR 
and predicts highly characteristic broad, double-peaked line 
profiles. In the present paper, we further develop this clumpy 
environment model and show that, taken with Paper I, there 
is reason to believe that these models can explain a signifi- 
cant fraction of UCHuR morphologies. 

In Section 2, we describe a model for UCHuR in which 
the dynamical effects of a central stellar wind are negligi- 
ble, so that the flow generated by mass injection is subsonic 
throughout the ionized region. This model is mainly appro- 
priate for early B and perhaps late O main-sequence stars 
which have less wind momentum than the massive OB su- 
pergiants assumed to be the exciting stars in Paper I. Mass 
loss from clumps is driven by photoionization and/or hy- 
drodynamic ablation. The former gives a mass injection rate 
that is independent of flow speed whereas, in the latter case, 



load ing is suppressed at low Mach numbers (Hartquist et al 



198f). We examine these two possibilities separately and al- 



low for possible radial variations in the mass loading rate 
by including a power-law dependence of the mass injection 
rate. Density and velocity plots, and line profiles and emis- 
sion measures are produced for several cases. We highlight 
the similarities and differences between the flows produced 
by the two different mass loading laws and show that the 
uniform spherical UCHuR morphology is naturally repro- 
duced along with the centre-brightened core-halo type. 

In Section 3 a hybrid model is described. A supersonic 
central stellar wind source is mass loaded to such an extent 
that it shocks and becomes subsonic before reaching the re- 
combination front. We produce an example of a line profile 
and emission measure plot for the case of a shock occur- 
ring close to the recombination front. We find that a shell 
morphology is produced if the emission from the high-speed, 
low-density supersonic gas is negligible. The velocity struc- 
ture in this case differs completely from that of the models of 
Paper I which have a similar overall morphology, thus pro- 
viding an observational test between the two. We note that 
it is likely that shell morphology UCHuR resulting from a 
fully supersonic mass loaded wind (Paper I) will in reality 
appear more clumpy than those that result from the partly 
subsonic structures described here. 



2 SUBSONIC AND TRANSONIC UCHuR 

We assume that the flow is subsonic throughout the ionized 
region. If the mass loading is dominated by photoionization, 
the rate is Mach number independent, while, if hydrody- 
namical ablation is the dominant mechanism, phenomeno- 
logical argume nts suggest that the vo lume mass loading rate 
varies as M 4//3 ( Hartquist et al. 1986). We discuss these two 



cases separately, but one should note that there may be more 
than one mode of mass injection at work in any given flow. 
The mass injection rate is also given a dependence on ra- 
dial distance in both cases. We ignore gravitational effects 
since 2GM,/c 2 <, OMR for an M» ~ 20 M star within an 
ionized region of characteristic radius R = 10 17 cm and an 
isothermal sound speed c ~ 10 kms _1 . We also assume that 
the flow is dust free, deferring a consideration of the dynam- 
ical effects of radiation pressure on a dust-gas mixture to a 
later paper (Williams, Dyson & Redman, in preparation). 
Finally, we assume that the dominant source of mass injec- 
tion is from the clouds. We take a volume mass injection 
rate 



q = q o M \r/R i ) c 



(1) 



where r is the radial coordinate, Rf is the distance of the 
recombination front from the star and a and f3 are constants. 

The continuity and momentum equations for isothermal 
flow are respectively 



dr 



(r pu) 



r q, 



(2) 



du c dp qu 
u— + --f=-2-, (3) 
dr p dr p 

where u and p are respectively the flow density and velocity. 
With the definitions 



M=-,r = ^— and P = P^--, 
c R f q Rt 



equations M and (H) g 



1 ve- 



il/ 



AM 



df (M 2 - 1) 



'M^- 1 {M 2 + 1) 



2M< 



dp_ 

df (M 2 ~ 1) 



(4) 

(5) 
(6) 



2.1 Mach number independent mass injection 

As noted above, this case is appropriate to the injection of 
material by photoionization, so that (3 = and equations Q 
and (^) give 



Mpr 2 = 



(3 + a)' 

For a ^ — 1, the velocity and density distributions are 



Ar = 



|(1 + a) + (3 + a)M 2 P 
1 



A 1+a |(1 + a) + (3 + a)M 2 |" (3 + a) 
where 



(7) 
(8) 

(9) 
(10) 
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Figure 1. The variation of the Mach number with fractional Figure 3. Emission measures as a function of fractional offset, 

offset z = r/Rf. From the lowest to uppermost curves, these cor- z. The a = —3/2 plot is the most centre-brightened, while the 

respond to a = 1, 1/2, 0, — 1/2, — 1 and — 3/2. a = 1 plot is the least centre brightened. 




0.4 0.6 
Offset z 



Figure 2. The variation of p, in units of c(3 + a)/qoRf, with z. 
From the lowest to uppermost curves (looking at the left hand 
side of the figure), these correspond to a = 1, 1/2, 0, — 1/2, — 
1 and - 3/2. 



2 + a 



and 8 - 



2 + a 



' (3 + q)(1 + q) ~ 3 + q' 

In the case a = — 1, equations ^ and (Jl^) become 

Af = Af- 1/2 exp(-l/(4M 2 )) and p = l/2M. (11) 

A is determined by the Mach number at the recombination 
front, i.e. by M = Mr. If Mr < 1 then the UCHuR is 
pressure confined by the external medium. If the external 
pressure is small, which we adopt as the most probable case, 
then the flow becomes transonic at the recombination front 
(cf. Williams & Dyson 1994) so that A = (4 + 2q)~ 7 (a / 1) 
and A = exp(-l/4) (a = -1). 

The Mach number and density distributions are given 
in Figs ^ and | respectively, as a function of the fractional 
offset z — r/R{, for six different values of a (including a — 0, 
which is uniform mass loading) . The densities are normalized 
such that the edge density is unity for each a (i.e. so p = 
1/(3 + a) at M = f — 1). 

The emission measure, EM = J n 2 dl, is given in Fig. ^ 
as a function of the fractional offset, z, from the central 
star, for each a. The individual plots are calculated using 
the density distributions of Fig. g. The a — —3/2 and —1 
distributions are centre-brightened and reproduce a core- 
halo morphology. A pho toevaporating disc, as described by 
Hollenbach et al. (1994), could provide additional mass in- 




Mach Number 



Figure 4. The profiles are almost identical but broaden slightly 
with decreasing a so that the outermost profile corresponds to 
a = —3/2. The dashed curve, shown for comparison, has no 
Dopplcr broadening. The fractional offset z = 0.1. 



jection close to the central star. However, on its own, a pho- 
toevaporating disc gives an approximately r~ 2 fall-off away 
from the disc, so additional mass injecting sources would be 
needed away from the centre. The other plots show progres- 
sively a more uniform brightness with increasing a. These 
other models could reproduce satisfactorily UCHuR that 
have spherical uniform morphologies. 

The corresponding optically thin recombination line 
profiles are shown in Fig. W. We assume that the gas emits 
with a Gaussian p rofile, and we neglect effe cts such as pres- 
sure broadening (Roelfsema & Goss 1992) which will add 

They 



characteristic Lorentzian wings to the profiles. They are 
viewed through lines of sight with offsets of z — 0.1 from the 
central star and are individually normalized to their peak in- 
tensity. The a = —3/2 is distinct but is only slightly broader 
than the higher a plots which are practically identical and 
appear almost superimposed in the figure. This shows that, 
in this case, line profiles will not be able to differentiate be- 
tween the models described above. These Gaussian-like pro- 
files are produced simply because the gas everywhere has a 
subsonic velocity. 
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Figure 5. Family of Mach number solutions plotted as a function 
of offset z for a = 0. We set M = 1 at the recombination front 
at z = 1 as one boundary condition and vary the inner boundary 
condition. 




Offset z 

Figure 7. The Mach number plotted as a function of offset z. The 
two extreme values of a are indicated. The initial Mach number 
is M = 0.1 at z = 0.2 and we also set M = 1 at z = 1. 
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Offset z 



Figure 6. The corresponding density solutions plotted as a func- 
tion of z for a = 0. The initial conditions are the same as in 
Figure 5. 



2.2 Mach number dependent mass injection 

We now take the mass injection rate to be proportional to 
M 4 / 3 which is appropriate for mass loading by hydrody- 
namical ablation. Substituting (3 = 4/3 into equations (El) 
and (h) gives 



dM _ / 
dr ~~ V 



M 



dp 
dr 



M 2 - 
2M 2 



2 M 1/3 (1 + Af 2 



{M 



/2 



M 



1/3; 



(12) 



(13) 



There is no non-trivial solution for hydrodynamic ablation 
that passes through f = 0, so we assume that the flow is 
started off by another mechanism such as photoionization, 
or by a stellar wind shocking, and match this at some radius 
to our model here. As one of our boundary conditions we 
have M = 1 at f = 1. Varying the other boundary condition 
(e.g. p at f ) generates a family of solutions for the velocity 
and density for each a. We present these solutions in Figs a 
and ^ for the case of a = only, for clarity. To show the 
variation of these solutions with a, we select the solution 
that passes through the point at f = 0.2 with M — 0.1 as 
an example. We retain the other boundary condition M = 1 
at r = 1 and examine the same six power laws as before. 
The full velocity and density solutions for this particular 





1 1 1 1 1 1 1 1 




a=-3/2 




'if 









0.4 0.6 

Offset z 



Figure 8. The density plotted as a function of z. The initial 
conditions are the same as in Figure 7. 



pair of boundary conditions are shown in Figs [7] and |^. The 
densities are normalized to the same edge density as before. 
The variation of the flows between different values of a is 
less marked than in the Mach number independent models 
of the previous subsection. The Mach number dependence of 
the mass injection rate means that mass loading is enhanced 
in a shell towards the edge of the region. This effect is most 
pronounced for the positive values of a. Gravitational ef- 
fects in fact may become important at the parts of the flow 
with particularly low velocities or at small radii, and we will 
discuss this elsewhere. 

In Fig. ^ we show the emission measure as a function of 
fractional offset z. Since there may be several ways in which 
the flow is initiated, we neglect any emission from this inner 
zone which is not part of the solution. This gives the dip for 
offsets less than z — 0.2, since our line of sight contains a 
sphere with no emission. It can be seen that the emission 
measures are quite similar for the different radial fall-offs. 
Comparing these with the emission measures for the flow 
loaded by photoionization ablation (shown in Fig. ^), we 
see that the Mach number dependent models are similar to 
the photoionization loaded flows with a <; 0, and thus they 
could also describe a UCH11R with uniform morphology. An 
optically thin line profile is shown in Fig. [l(^ for a fractional 
offset z — 0.2. The profile is hard to distinguish from the 
profiles of the photoionization loaded flows with a ^ seen 
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Offset z 

Figure 9. Emission measures for the same initial conditions and 
radial dependences as above. The radius at which the flow begins 
is indicated by the dashed line. 
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Figure 11. Emission measure as a function of z for a supersonic- 
subsonic flow. The position of the shock is indicated by the dashed 
line. 




Mach Number 



Figure 10. An example of a line profile for the same initial con- 
ditions and radial dependences as above with z = 0.2. The six 
profiles are practically identical so we show just one here. The 
dashed curve, shown for comparison, has no Doppler broadening. 



in Fig. ^| Only one profile is shown because the six models 
give practically the same line profiles. 



3 UCHnR IN TERMS OF A 

SUPERSONIC-SUBSONIC FLOW WITH AN 
INTERNAL SHOCK 

In this section, we combine the results of the previous sec- 
tions with those of Paper I to allow for the possibility that 
the flow is decelerated by mass loading to such an extent 
that a wind termination shock occurs before the recombi- 
nation front is reached. This will occur if the Mach number 
in the supersonic flow is predicted by the ballistic approxi- 
matio n to reach values of M 2 before th e recombination 
front (Williams, Hartquist fc Dyson 1995). In Paper I we 



used the conservation of mass and momentum to find veloc- 
ity and density distributions for a fully supersonic, uniformly 
mass loaded flow: 



M 



3/i* 



P = 



Anq 2 r A 



(14) 



The mass loading mechanism in the supersonic fl ow before 
the shock r adius is independent of Mach number ( Hartquist 
et al. 1986). The emission will be dominated by the low- 



velocity, high-density gas at the edge of the mass loading re- 
gion. The high velocity of the unshocked wind will mean that 
it stands well clear of the near-Gaussian profiles predicted 
for the subsonic region. However in all but the cases with 
the lowest pre-shock velocities, or thinnest subsonic shells, 
the high density in the subsonic region (which is a factor 
substantially more than 4 higher than the pre-shock density 
because of the cooling of the shocked gas) will mean that 
this low-velocity gas dominates the emission. For simplicity, 
we assume that the mass injection rate is independent of 
flow velocity and radial distance so that /3 = a = 0. 

In the emission measure plot shown in Fig. |ll| a shell 
morphology is clearly apparent. The emission at offsets less 
than the shock radius is from the subsonic gas in our line 
of sight and not from the supersonic gas which contributes 
very little to the emission. In fact, even if there were notice- 
able emission from the supersonic gas, the overall emission 
measure would still be essentially that of a shell because 
the density distribution (equation |l4| ) is strongly weighted 
to the edge of the supersonic region. In Fig. h2j we show 
an optically thin line profile for this shock radius taken at 
a line of sight with offset z — 0.1, close to the centre. The 
profile consists of two near-Gaussian components, generated 
by the parts of the shell that are moving towards and away 
from us, separated in velocity by an amount that is a frac- 
tion of their widths. This results in a line profile that is 
broader than those produced using the fully subsonic mod- 
els described earlier. We find very similar results if the mass 
loading is Mach number dependent. 



4 CONCLUSIONS 

In this paper we have modelled UCHnR as subsonic flows 
from photoionized clumpy clouds lying in the vicinity of an 
ionizing stellar source. We have found that, by assuming a 
simple spatial distribution of clumps and allowing the flow 
to be supersonic, subsonic or both (with a separating shock), 
we can reproduce three different UCHnR morphologies. 

For the fully subsonic models, we have found that a 
radial dependence of the mass injection rate is more im- 
portant for the photoionization-induced subsonic flow. The 
emission measures show that the M 4 ^ 3 flows are less centre- 
brightened than the photoionization flow. Either of the two 
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dust will have on these models, consider the intermediate- 
scale structure in the regions and incorporate any gravita- 
tional effects. 
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Figure 12. Line profile for z = 0.1 taken through a region with 
a shock at 0.9-Rf. The dashed curve, shown for comparison, has 
no Doppler broadening. 



types of model can naturally describe the uniform spher- 
ical UCHllR. In addition, the photoionization flow with 
q = —3/2 or —1 can satisfactorily explain the core-halo 
UCHllR. Line profiles are of no use in distinguishing be- 
tween the two different flows. For both cases we have calcu- 
lated velocity and density distributions that will be used in 
more detailed future work. 

Paper I showed that a supersonic wind-blown model 
could reproduce the shell morphology but with broad, 
double-peaked line profiles. In contrast, we have suggested 
in Section 3 that if a shock forms before the recombination 
front then a shell morphology is retained but with a sin- 
gle peaked line profile. Which type an individual object is 
will be settled by high-resolution observations of the velocity 
structure of the region. 

At high enough resolution, the largest individual clumps 
will begin to be resolved. The irregular/multiple-peaked 
UCHllR are perhaps simply a collection of large clumps 
which are ionized on their outsides by the central source. 
The 'hyperco mpact ' continuum sources found in clusters by 
Gaume et al. ( 1995 ) are perhaps clumps such as this. Prelim- 
inary work suggests that a non-spherically symmetric distri- 
bution in the mass loading around the region can describe 
the arc- like UCHllR (Paper III). 

The solutions to the subsonic flow equations may be 
applicable to other objects. Although visible H II regions are 
older and more evolved than UCHllR, their morphol ogies 



hav e be en found to be very similar to those of UCHllR ( Fich 
1993|)71md this may imply that they could be used to test 
some of the observational diagnostics described here and in 
Paper I. The extent of the neutral H I envelopes around H n 
regions is explained by t he structure of pho todissociation re- 
gions in clu mpy cl ouds (Howe et al. 1991). However, Roger 
& Dewdney (1992) note that the outflowing Hi shocks seen 
around Hn regions (e.g. Kuchar & Bania 1993) are not as 
thin and do not have as high densities as predicted by stan- 
dard theory. Such low-density shells of Hi may result from 
support by the swept-up magnetic fields ( [Mathews 



O'Dell 



196S()~3Ve would suggest that the neutral winds beyond re- 
combination fronts would also explain these structures. If 
this is the case, then observations of Hn regions agreeing 
with our predictions for UCHllR would lend support to our 
models. 

We will address elsewhere the dynamical effects that 
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